Electronic and optical properties of Cadmium fluoride: the role of many-body effects 
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Electronic excitations and optical spectra of CdFo are calculated up to ultraviolet employing state-of-the-art 
techniques based on density functional theory and many-body perturbation theory. The GW scheme proposed 
by Hedin has been used for the electronic self-energy to calculate single-particle excitation properties as en- 
ergy bands and densities of states. For optical properties many-body effects, treated within the Bethe-Salpeter 
equation framework, turn out to be crucial. A bound exciton located about 1 eV below the quasiparticle gap is 
predicted. Within the present scheme the optical absorption spectra and other optical functions show an excel- 
lent agreement with experimental data. Moreover, we tested different schemes to obtain the best agreement with 
experimental data. Among the several schemes, we suggest a self-consistent quasiparticle energy scheme. 
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I. INTRODUCTION 

Fluorides and fluorite-type crystals have recently attracted 
much interest for their intrinsic optical properties and their 
potential applications in optoelectronic devices. Because of 
their peculiar optical properties, possible applications in ul- 
traviolet (UV) laser optics could be considered. For example, 
excimer ArF lasers with an emission wavelength of 193 nm 
and KrF lasers with an emission wavelength of 248 nm offer 
many applications in medicine: e.g., photoablation of the sub- 
stantia propria (or stroma of cornea), and high precision tissue 
ablation. This comes in addition to the common application 
in photolithography of the semiconductor industry. 

Here we will consider cadmium fluoride (CdF^) which is 
characterized by a large gap energy (of the order of 9 eV) 1 
and, hence, a high transparency in a wide energy range, while 
the more famous calcium fluoride (CaF-i) has a direct band 
gap at r of 12. 1 eV and an indirect band gap estimated around 
11.8 eV. 1 CdFi has been chosen as a representative candidate 
for a material class which seems to exhibit similar behavior, 
as well as numerical problems (e.g., BaF2 ). In particular, its 
overall very small values of the dielectric function, the ex- 
tended high-energy tail, and the full ionicity (i.e., the weak 
band dispersions) may pose critical theoretical and numerical 
problems. 

By doping, usually with trivalent metal impurities and a cer- 
tain thermochemical treatment, CdFj can be also produced in 
the semiconducting state. Many dopants, as Sc and Y, gen- 
erate only a shallow donor state, but In or Ga related de- 
fects exhibit a bistable behavior. In addition to the shallow 
state these impurity centers also possess a deep state.' 1 " 8 Ap- 
plication of UV and visible light at low temperatures results 
in disappearance of the absorption peak corresponding to the 
deep electronic state and an infrared absorption band associ- 
ated with the occupation of the metastable shallow state. In 
other words, the electrons occupying the deep state are lifted 
due to photoexcitation to the shallow state. At low tempera- 
tures the electrons cannot be trapped back to the deep states 



due to a barrier induced by different atomic configurations for 
shallow and deep states. The change in the electronic oc- 
cupation causes a large difference in the local refractive in- 
dex, which can be used as means for holographic writing at 
nanoscale spatial resolution. A general observation is that 
the formation energies of most defects are found to be very 
low. This suggests that the fluorite structure in CdFi often 
contains high defect concentrations in agreement to what is 
found experimentally/ It is also essential that CdFi crystals 
of quite large size can be obtained at a relatively low cost. 

Optical properties of CdFj have been experimentally deter- 
mined by different spectroscopic techniques since the seven- 
ties of the previous century 11 ' 1 ' in the fundamental absorp- 
tion region and in the core-level excitation range. In the same 
decade reflection spectra of CdF2 have been determined in 
comparison with those of SrFi, 1 while UPS and XPS spectra 
of CdFi and SrFi have appeared in the literature in 1980. 
Also j3-PbF2 and CdF2 mixed crystals absorption coefficients 
have been reported by spectrophotometry measurements. 14 

Various theoretical methods have been applied to study ei- 
ther the ground state or the excited states of the fluorite com- 
pounds. The energy bands and reflectance spectra of CaFi and 
CdFi have been determined within a combined tight-binding 
and pseudopotential method. Mixed crystals of CaF%, SrF%, 
CdFi, ji-PbFz have been studied with respect to their elec- 
tronic energy bands and density of states (DOS) within the 
linear muffin-in orbital (LMTO) method. Linear and non- 
linear optical properties of the cubic insulators CC1F2, SrF2, 
CdF2, BaF% and other compounds have been determined by 
first-principles orthogonalized linear combination of atomic 
orbitals (OLCAO). Point defect studies in CdFi have been 
performed within the plane wave pseudopotential (PW-PP) 
method. With respect to one-particle and two-particle elec- 
tronic properties and energy band gaps, state-of-the-art tech- 
niques have been applied until now only to CaFi. In feet, elec- 
tronic band structures of CaFi have been determined within a 
GW approximation, using a PW-PP scheme. On the other 
hand the imaginary part Im e(o) of the dielectric function has 
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been calculated for CaF-i after an iterative procedure using an 
effective Hamiltonian, within a PW-PP scheme considering 
a screened interaction for electron-hole (e-h) coupling, and 
using localized orbitals. 

In this paper the electronic excitations and optical spec- 
tra of CdFi are calculated up to 40 eV employing state-of- 
the-art techniques based on density functional theory (DFT) 
and many-body perturbation theory (MBPT). We use the GW 
scheme proposed by Hedin -1 for the exchange-correlation 
(XC) self-energy, employing an Heyd-Scuseria-Ernzerhof 
(HSE) hybrid functional as starting point for the electronic 
quasiparticle structure to calculate single-particle excitation 
properties as the energy bands and the DOS. The DOS near the 
gap region and the energy-band structure are compared with 
existing data from literature; we show and discuss the com- 
parison with them. The role of many-body effects turns out to 
be fundamental for these single-particle properties. Moreover, 
we demonstrate that a simple one-shot GW (Go Wo) is not suf- 
ficient to obtain the correct band gaps. For optical properties 
excitonic effects, treated within the Bethe-Salpeter equation 
(BSE) framework, are crucial to allow a reasonable com- 
parison with existing experimental spectra as well. We dis- 
cuss the electronic excitation structure and the optical spec- 
tra including the existence of bound excitons, and compare 
with available experiments. One of our goals is to provide a 
wide scenario of the one- and two-particle properties of this 
material treated within first-principles techniques. However, 
as considered later on, a very extended high-energy tail and 
overall rather small values of the dielectric function require 
inclusion of an unusually large number of bands or band pairs 
in the GW and BSE schemes. In regard to this last point, 
we suggest a theoretical framework and a simulative protocol 
able to give solutions with appreciable agreements with the 
available experimental measures, saving at the same time the 
computational demands. 



II. GROUND STATE PROPERTIES AND 
COMPUTATIONAL DETAILS 

All the calculations have been performed using density 
functional theory (DFT) - as implemented in the plane- 
wave basis code VASP. ' The projector augmented wave 
(PAW) ' method is applied to generate pseudopotentials 
and wave functions in the spheres around the cores. In stan- 
dard use, VASP performs a fully relativistic calculation for the 
core-electrons and treats valence electrons in a scalar relativis- 
tic approximation. 25 ' 26 ' 29,30 The effect of spin-orbit coupling 
on the energy bands has been considered extensively in Ref. 
31, where the spin-orbit effect, in the case of C1JF2, slightly 
decreases the gaps by about 0.05 eV only. Thus spin-orbit 
effects will be neglected in the present paper. 

We performed the calculations using different XC func- 
tional. The local density approximation (LDA) has been 
used for the XC energy, as given by Ceperley and Alder, 
and parametrized by Perdew and Zunger. • In addi- 
tion, the generalized-gradient approximation (GGA) in the 
parametrization of Perdew, Burke, and Ernzerhof (PBE) is 



used. 14 Moreover, we considered also a revised version of the 
PBE functional which gives better results in solids which we 
shall refer to as PBEsol 35 from now on. 

The lattice structure of CdF^, as well as of all the other 
fluorides, is a cubic one with the space group Fniim, with 
three ions per unit cell, i.e., one cation Cd placed in the origin 
and two anions F situated at ±(|a, |a)/ 6 In the crystal 
the F^-' ions form a simple cubic sublattice surrounded by 
a face-centered cubic lattice of Cd( ++ > cations. All fluorides 
with cations belonging to the II and IIB groups are stable in 
this crystallographic structure. 

Besides the Cd(5s) and F(2s,2p) valence states also the 
shallow Cd(4s, 4p, Ad) core states are treated as valence elec- 
trons. Inclusion of the Cd(4d) states is important because 
they appear in between the F(2s) and F(2p) bands. Inclu- 
sion of the strongly bound Cd(4s,4p) states (about 100 or 60 
eV below the valence band maximum) is not an absolute must 
but helps a lot to construct s- and /9-pseudopotentials with ex- 
cellent scattering properties up to very high energies which 
are needed for a proper treatment of a large number of un- 
occupied bands in the calculation of the dielectric function. 
The wave functions are expanded in a plane-wave basis set 
with a cutoff energy of 950 eV. The DFT/HSE calculations 
for this cutoff value are fully converged. For GW/BSE calcu- 
lations the error bar has been estimated at about 0. 1 e V. The 
face-centered cubic (fee) Brillouin zone (BZ) is sampled by 
F-point centered meshes with 16x16x16 Monkhorst-Pack k- 
points (converging within a few meV). The minimum of the 
total energy with respect to the volume is obtained by fitting 
to the Vinet equation of state. For the calculation of the co- 
hesive energies £«,/,, we have subtracted the spin-polarized 
ground-state energies of the free atoms. Comparing exper- 
imental measurements with results calculated with differ- 
ent XC functionals, see Table I, the best results are obtained 
within the PBEsol scheme. Therefore we have decided to use 
from now on the PBEsol lattice constant of 5.39 A consis- 
tently for all calculations (ground-state, one- and two-particles 
excitations). This is an important point since band structure 
test calculations on DFT level for different lattice constants 
show that the gap value changes by about +0.1 eV (-0.1 eV) 
when going from the PBEsol to the LDA (PBE) lattice con- 
stant. 



TABLE I: Structural data for CdF~i are given: lattice constant 
a Q , bulk modulus B and its pressure derivative B' Q = dB Q /dP, 

and cohesive energy E co f,. Three local and semilocal XC 
functionals are used, as discussed in the text. Other 
theoretical results and experimental references are reported in 

column "OTheory" and in column "Expt.", respectively. 





LDA PBE 


PBEsol OTheory 


Expt. 


flo (A) 


5.30 5.49 


5.39 5.39 lh 


5.36-5.39 


Bo (GPa) 


126.70 93.80 


108.50 123.00 16 


114.60 39 


B' a 


4.76 4.88 


4.85 4.85 




Ecoh (eV) 


-12.07 -9.55 


-10.39 
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m. QUASI-PARTICLE EXCITATIONS 

The Kohn-Sham (KS) eigenvalues of the DFT scheme can- 
not be interpreted as energies of single-particle/quasiparticle 
(QP) electronic excitations and, therefore, cannot be com- 
pared with the band structure and the DOS from experi- 
ments. We apply the many-body perturbation theory within 
Hedin's GW approximation for the electronic self-energy 
operator." 1 It can be demonstrated that one has to solve the 
QP equation: 23 - 40 " 42 

{-^V? + v fl (r) + v«(r)}*f (r) 

+ J «iVl(ry; £ f)<(r') = ef *f (r). (1) 

Equation (1) has a similar structure as the KS equation of DFT. 
The fundamental difference is that the XC potential of DFT 
is replaced by the non-Hermitian, nonlocal energy dependent 
self-energy operator E resulting in QP eigenvalues e9 F and 
QP wave functions < t , ^ P - Due to the fact that the Hamiltonian 
is non-Hermitian, QP eigenvalues are usually complex num- 
bers. The real part describes resonance energies (excitation 
energies) and the imaginary part defines the lifetime of ex- 
cited states. Equation (1) should be solved self-consistently; 
however, the most common GW approach avoids this proce- 
dure. Since empirically one often finds that QP wave func- 
tions are similar to the KS DFT-LDA ones, it is natural 
to take them as a starting point for a GW calculation and, 
hence, to calculate QP corrections in the sense of a perturba- 
tion treated within first-order perturbation theory. 2 ' 40 ' 41 The 
following equation, 

ef = £k + (<P*|£(ef ) - VxcWk) , (2) 

is used in conjunction with the random phase approximation 
(RPA) dielectric function to calculate the screened Coulomb 
interaction Wo and the DFT Green's function Go- This ap- 
proach called GoW) usually works very well and gives very 
good results for the GW eigenvalues (band gaps, bandwidths 
and band dispersions). 40 ' 41 Another GW treatment of the 
QP band structure is based on the so-called generalized KS 
schemes (gKS). The gKS methods have in common that 
nonlocal exchange-correlation potentials involving partial or 
screened exact exchange are inserted into the KS equations 
combined with (semi-)local DFT potentials. From the view- 
point of band structure calculations these gKS schemes can 
be also considered as oversimplified GW schemes. Hence, 
they can be considered as an improved starting point (eigen- 
values, wave functions) for Go Wo calculations. While results 
are similar for many systems, gKS schemes as a GoW) start- 
ing point can provide much improved results for those cases 
where DFT-LDA or -GGA fails drastically (wrong band or- 
dering, extreme gap underestimates). Since the experimental 
gap is almost three times as large as the DFT-PBEsol one (a 
gap correction of about 5 eV is necessary to arrive at the ex- 
perimental values) a gKS starting point (still leaving a 3 eV 
gap error) can be considered as an improvement also for fluo- 
rides. In the present work QP bands of CdF2 are determined 
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FIG. 1 : (Color online) Energy bands of CdF2 calculated 
within the Go Wo (light lines) and the scQP-GW (dark lines) 

schemes. On the right panel of the figure, the DOS are 
reported for valence bands (VB) and conduction bands (CB). 
Notice the magnification factor used for the conduction 
bands. 



by an iterative solution of the QP equations given in Eq.(2) 
with the XC self-energy in the GW approximation. The itera- 
tion starts with the gKS equations, with a self-energy derived 
from the nonlocal HSE03 hybrid functional (zeroth order). 22 
Subsequently, the GW corrections are calculated in first-order 
perturbation theory, i.e., within the one-shot GoWo. Thereby, 
the full frequency dependence of the RPA dielectric function 
entering the screened Coulomb interaction (Wo) is taken into 
account sampled on a frequency grid with 128 points up to 
H(0 «350 eV, i.e., no plasmon pole models 42 or model di- 
electric functions 1 1 have been involved. We present here QP 
eigenvalues calculated in the GoW) approach on top of the 
HSE03 ground-state electronic structures. 

The strategy to calculate QP corrections then follows: The 
band structure and DOS have been calculated first with the 
PBEsol GGA, then a HSE03 calculation is performed, and fi- 
nally we apply the GoW) scheme on top of HSE03. Even after 
this step, the QP energies are of the order of 1 eV smaller than 
the experimental ones. For comparison we made also Hartree- 
Fock (HF) calculations, which show as expected deviations in 
the opposite direction (i.e., an unrealistic overestimate of the 
gaps). • As the last step, in order to open that gap further, 
we performed a "self-consistent quasiparticle energy" calcula- 
tion within the GW approximation (namely, scQP-GW) which 
is a kind of "iterated GoW)" just updating the eigenvalues only 
but keeping the HSE03 wave functions fixed. This procedure 
gives almost perfect gaps with respect to experiment. ' 47 

The main outcomes of the present paragraph are reported 
in Table II and in Figs. 1 and 2. More details are also dis- 
cussed in Ref. 48. The QP values for the fundamental gaps 
F F (direct) and W — F (indirect) obtained within different 
approximations are also reported in Table II. It is clear from 
Table II and from Figs. 1 and 2 that after the computational 
procedure adopted here an excellent agreement has been ob- 
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FIG. 2: (Color online) Density of states around the gap 
region calculated within different approximations: HF, 
PBEsol, HSE03, G W Q , and scQP-GW. Notice the 
magnification factor (five) used for the conduction bands. 



estimate resembles a similar behavior taking place in other 
systems. ' Notice that the ^-electrons of Cd play an impor- 
tant role. Indeed, they are present inside the valence bands 
in the energy range displayed in Fig. 1, and give rise to flat 
bands/sharp peaks near -5 eV, due to the localized character 
of these states. Localization of CB states is much weaker, and 
hence their dispersion is much stronger. Therefore, the aver- 
age height of the DOS is smaller in comparison with the VB 
states. It is interesting to report in detail the energy band gaps 
calculated within the two different GW schemes proposed in 
the present publication, namely Go Wo and scQP-GW. As it is 
clear from Table III a nearly rigid shift of about 1 eV appears 
going from GoW) to scQP-GW transition energies. Therefore, 
with respect to CdF%, considering Fig. 1 and Table III, the 
action of the self-consistency upon the eigenvalues operates 
a kind of "second order scissor operator rule" on the conduc- 
tion bands to obtain good accordance with experiments. The 
valence bandwidths after the two approximations of XC, re- 
spectively GoW) and scQP-GW, agree within 0.5 eV. Thereby, 
the scQP-GW value of 22.99 eV corresponds excellently to an 
experimental estimate of 23 eV reported by Raisin. 3 



tained with available experiments. 13 ' 49 The other point is that 
for CdF2, Go Wo is not sufficient and a self-consistent proce- 
dure, GW based, has to be considered in addition. Referring 
to Fig. 1 the main discrepancies between the different GW 
methods appear in the conduction bands. This point is also 
clear from the DOSs which are nearly overlapping for the VB 
energy region, as shown in Figs. 1 and 2. 

Therefore, the electronic structure of CdFi, within the 
PBEsol, Go Wo, and scQP-GW schemes, shows an indirect gap 
at W — r of 7.83 eV (expt. 7.8 eV 4 '') and a larger direct gap 
at r of 8.58 eV (expt. 8.4 - 8.7 eV). 13 - 49 The best accordance 
with experiments is found within the scQP-GW scheme with 
a high level of accuracy. In Fig. 2 we present the calculated 
DOS within different schemes. As is clear from this figure the 
VB regions after the HSE03 and PBEsol methods seem very 
similar while discrepancies arise for the CB region between 
the results after the two methods (higher value of the gap after 
the HSE03 scheme). 

Similarities for the VB regions appear also between the re- 
sult after the Go Wo and scQP-GW schemes. The latter method 
does operate a further opening of the gap. In the upper part of 
Fig. 2, the HF DOS is reported. It determines a (too) large 
overestimate of the band gaps (about 16 eV at F) and valence 
bands are only partially (the shallow ones) similar to those 
calculated after the other methods. This HF band gap over- 



TABLE II: QP energies for the fundamental gaps of CdF2, 
between lowest CB and highest VB at high symmetry points, 
are computed with several approximations and compared 
with available experiments. 



data in eV HF PBEsol HSE03 GqWq scQP-GW Exp. 
r-T 16.02 3AI 5J50 1M 8^58 8.4 ,8.7 
W-T 15.30 2.91 4.91 6.90 7.83 7.8 



IV. TWO-PARTICLES PROPERTIES AND DIELECTRIC 
FUNCTION 

The treatment of excitons and resulting optical spectra, i.e., 
the full excitonic problem, requires to set up and diagonalize 
an electron-hole Hamiltonian H. Within Hedin's GW scheme 
and a restriction to static screening this two-particle Hamilto- 
nian reads for singlet excitations in matrix form as.' 

6(vck,v'c'k') = [£f(k) - ^(k)]^^ 

-J d\ j d 3 r / ( P ; k (r)(p 6 , k ,(r)W(r,r')(Pv k (r')(p i * k ,(r') 

+2 J d\ J fl f 3 rV* k (r)(p vk (r)v(r,r')(p,'^(r / )(p; k ,(r / ), (3) 

where matrix elements between KS or gKS wave functions 
of CB states and VB states occur. Contributions to Eq.(3) 
which destroy particle number conservation have been omit- 
ted. The first term describes the noninteracting quasi-electron 

TABLE III: Transition energies between lowest CB and 
highest VB at high symmetric points, and valence bandwidth 
within the two GW schemes. All data are expressed in eV. 



Direct bandgaps 


G W 


scQP-GW 


L- 


■¥ L 


12.11 


13.08 


r- 




7.64 


8.58 


X - 


*x 


12.57 


13.62 


w - 


+ w 


12.28 


13.66 


K- 




12.65 


13.44 


Indirect bandgaps 




W - 




6.90 


7.83 


W - 




12.50 


13.55 


W - 




11.56 


12.54 


W - 




12.47 


13.44 


VB width (F 2s - F 2p) 


22.47 


22.99 
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quasi-hole pairs. The second term accounts for the screened 
electron-hole Coulomb attraction with the statically screened 
Coulomb potential W(r,r'). The third contribution, gov- 
erned by the nonsingular part of the bare Coulomb interac- 
tion v(r,r'), represents the electron-hole exchange or crystal 
local-field effects. ' The QP eigenvalues entering the first 
term and W(r,r') entering the second term have been directly 
extracted from the scQP-GW results. After diagonalization 
of the exciton matrix, or, more precisely, solving the homo- 
geneous BSE or stationary two-particle Schrodinger equation 
with the eigenvalues E A and the eigenfunctions A A (vck) of 
the pair states A, frequency dependent macroscopic dielectric 
function including excitonic effects can be written as: 



Energy (eV) 
15 20 25 



167te 2 fi 2 ^ 
£aa((0) = O aa H ~ 2- 



V 



1 



L(ck|v a |vk) 
„ k iw^i) Mvck) 

1 



E A -h(a> + iY) E A +h(co + iy) 



(4) 



where \> a is the corresponding Cartesian component of the 
single-particle velocity operator and y is the pair damping 
constant. The crystal volume is given by V. The details of 
the standard scheme using the direct diagonalization of Eq.(4) 
have been discussed elsewhere. 18 ' 3 ' 50 ' 5 

Since the rank of the Hamilton matrix in Eq.(4) is extremely 
large (it is given by the number of valence bands times the 
number of conduction bands times the number of k points) 
a straightforward diagonalization of this matrix is often not 
possible due to high CPU but also memory requirements. 

Therefore, we mainly refer to a numerically efficient 
scheme 4 to solve the Bethe-Salpeter equation, which has 
been successfully applied by Bechstedt and coworkers on sev- 
eral systems. This method is based on the calculation of the 
time evolution of the exciton-state by Fourier transform on 
the time domain. This scheme delivers optical spectra directly 
but no exciton eigenvalues and eigenfunctions. hl The calcula- 
tion of the co-dependent polarizability can be considered as an 
initial-value problem. This means that it can be written as 
a Fourier representation within an integral of time-dependent 
elements. The time evolution of these elements is driven by 
the pair Hamiltonian [see Eq.(3)] and it is operatively per- 
formed by using a central-difference method, which requires 
one matrix-vector multiplication per time step. The time in- 
tegration can be truncated due to exponential decay factors 
and it determines a scheme nearly independent from the di- 
mension N (number of pair states) of the system. Two fur- 
ther advantages come from the matrix-vector multiplication 
scheme, i.e., an 0(N 2 ) dependence on operations count and 
the possibility to effectively distribute the multiplications on 
several processors of a parallel computer. 46 Applications of 
this method to the calculation of the co-dependent dielectric 
function of bulk and surface systems have been successfully 
performed with good comparison either with experimental re- 
sults and with outcomes of the matrix diagonalization scheme. 
Moreover, a further speed-up of the calculations could be ob- 
tained combining this method with the use of a model dielec- 
tric function, ' ■ which, however, has not been employed in 
the current calculations which involve no models at all. This 




FIG. 3: (Color online) Imaginary and real part of the 
dielectric function calculated within the BSE scheme in 
comparison with experimental results after Ref. 12. 



scheme has the big advantage that it never requires the ex- 
plicit diagonalization of the Hamilton matrix, because it needs 
only to evaluate the Hamilton exciton matrix times the current 
wave function for every time-step. Massive CPU-time saving 
and improved scaling of the computational work with respect 
to system size result. On the other hand, we should admit that 
the lack of eigenvalues and eigenfunctions limits the detailed 
knowledge of excitons, e.g., their real-space representation. 
The investigation of the fine structure of optical absorption 
spectra, especially near the absorption edge and for low en- 
ergy optical transitions, requires finer fc-point samplings than 
those necessary for the ground-state calculations. Due to the 
large number of bands and high cutoffs involved we have, 
however, to restrict ourselves to a rather limited number of 
^-points. For the GW and BSE calculations, we apply a F- 
centered 8x8x8 mesh. Furthermore, a sufficiently large num- 
ber of conduction bands has to be included in the calculations 
to describe the electron-hole pair interaction properly. The 
number of empty states is limited by introducing a cutoff en- 
ergy for the electron-hole single-particle transitions (without 
QP corrections) that contribute to the excitonic Hamiltonian 
of Eq.(3). This transition energy cutoff was chosen to be 60 
eV, corresponding to up to 50 unoccupied bands. A rank of the 
excitonic Hamiltonian matrix of the order of 210.000 results 
from this setup. In addition, for the calculation of the full fre- 
quency dependent dielectric function, screened exchange in- 
tegrals, and four-orbit integrals in the GW and BSE schemes, 
one may introduce a lower plane-wave cutoff than for ground- 
state calculations without spoiling accuracy too much. In our 
case we use a cutoff of 475 eV. 

The matrix elements of the excitonic Hamiltonian have 
been evaluated by means of HSE03 wave functions. Since 
the HSE03+G()Wo QP eigenvalues still leads to a gap under- 
estimate (see Table II), HSE03+scQPGW eigenvalues are em- 
ployed instead of HSE03+GoWo eigenvalues on the main di- 
agonal of the Hamiltonian. It has been carefully checked that 
no change of band ordering occurs when going from PBEsol 
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Energy (eV) 

FIG. 4: (Color online) Imaginary part of the macroscopic 
dielectric function calculated within different 
approximations, namely PBEsol, HSE03, GqWq, scQP-GW, 
and BSE. Notice that each graph has been shifted by 20 units 
along the y axis. 



to HSE03. 

In Fig. 3, we report the principal result of the present calcu- 
lations: the imaginary part of the dielectric function Im e(ft)) 
determined within the BSE approach by using the method 
of Schmidt et al. . " The real part has been obtained by a 
Kramers-Kronig (KK) transform afterwards. Also the exper- 
imental results after Bourdillon - are reported, which are ob- 
tained with synchrotron radiation technique. The agreement 
both in the positions and intensities of the peaks between the- 
ory and experiment should be stressed for the imaginary part 
of the dielectric function (e.g., the first main three peaks are 
located at 7.7, 13.1, and 15.1 eV, respectively). As shown in 
the top panel of Fig. 3, the good agreement turns out for the 
full energy spectrum from near UV up to 40 eV, either below 
and above the QP direct gap. The first structure is a F exci- 
ton, while the second and the third peaks are excitons which 
may originate from X, L, or W transitions between valence and 
conduction bands. Since we calculated the spectra directly by 
the method of Schmidt, we have no access to excitonic eigen- 
values or eigenfunctions and can just speculate on the basis of 
single-QP eigenvalues. 

However, on the basis of the experimental spectra of 
Bourdillon as shown in Fig. 3, the value of &*, is seem- 
ingly about 1.8, which is in large discrepancy to our value of 
2.5. However, one should mention that the data of Bourdil- 
lon have been extracted indirectly from reflection spectra and 
that there exist other values for £«, from more direct meth- 
ods. Most probably, the best directly measured value is about 
2.4, and it is also indirectly confirmed by Krukowska x ^ from 
measurements of the refraction index at visible optical fre- 
quencies, obtaining a value slightly larger than 2.4, i.e., about 
2.46-2.5. Furthermore, by using the Lyddane-Sachs-Teller 
(LST) relation M ' the consistency of our result could also be 
indirectly confirmed from measurements of Sq (which yield 
values of 7.7, 7.8, and 8.3-8.5 at K, 80 K, and 300 K, 
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FIG. 5: (Color online) Real and imaginary part of the 
refractive index (n,k) as a function of energy in comparison 
with measured spectra (after Ref. 12). 

respectively) together with phonon data, 39 ' 58 which lead to- 
ward values in the range of 2.3-2.5 for £„. Hence all these 
values from direct measurements are in excellent agreement 
with our computed value. Furthermore the real part of the di- 
electric function Re e(o), as shown in the second panel of 
Fig. 3, is in close agreement with the experimental data after 
Ref. 12. 

In Fig. 4, we show the imaginary part of the dielectric func- 
tion Im e((o) within the different approximations proposed 
here, starting from the DFT reference which is the PBEsol 
one. Going from this approximation to the gKS HSE03 one 
an opening of the gap results. Both GW schemes further en- 
large the gap energies. Finally the BSE results are disposed 
which reproduce the experiment at best in terms of oscillator 
strengths, peak positions, and onset behavior. 

V. OPTICAL AND LOSS FUNCTIONS 

In Fig. 5, the imaginary and real part of the refractive index 
in comparison with experiment after Ref. are shown. A fair 
agreement results with respect to the three main peaks at 7.8, 
13.2, and 15.2 eV and to the high-energy region 20-40 eV. 

The reflectivity function R is show in Fig. 6 in compar- 
ison with an experimental spectrum . The positions of the 
main peaks of the theoretical curve (namely, occurring at 7.8 
eV, 13.2 eV, and 15.2 eV) reproduce the corresponding ones 
of the measurement, except for the fourth peak (at 18.0 eV) 
which is overestimated in intensity and it should be compared 
with the shoulder around 17 eV as obtained in the experimen- 
tal spectrum. Furthermore, the energy region above 16 eV 
is only in qualitative agreement with experiment. Regarding 
the first three peak positions of the reflectivity R, our results 
compare well with the experimental values after Berger and 
coworkers. 1 ' Albert et al. 15 assign these peaks to rX — > Tj c , 
XZ — > X^ and XY, —> X^ excitons. The first excitonic peak 
at 7.7 eV is consistent with the results of 7.6 eV given by 
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FIG. 6: Calculated reflectance spectrum within BSE and 
experiments after Ref . 1 2 as a function of energy. 



1.2 I 1 1 1 1 1 1 

[ BSE 

Exp. 

1.0 - 




5 10 15 20 25 30 35 40 



Energy (eV) 

FIG. 7: Electronic energy loss function after present 
calculation vs experiments (Ref. 3). 



Orlowski 4 ', Raisin' 3 , and by Forman and coworkers. 

We focus on the fact that in our computations a natu- 
ral/instrumental broadening of 0.2 eV has been applied in all 
the curves. Within our scheme, "our experimental curve" is 
the imaginary part of the dielectric function, Im e(co) (the first 
outcome of the computational code). All the other functions 
are derived from this calculation. Thus, an important method- 
ological point is the following: From each optical function 
produced here one can return to the original dielectric function 
which was the output of the code we used ("our reference re- 
sult"). Moreover the comparison with the experimental result 
for Im e(ffl) is anyway excellent (see Fig. 3). If one operates 
a broadening to fit at best an optical function with respect to a 
particular experiment, the reverse operation (to generate again 
the Im e(o) in agreement with experiment) is not guaranteed. 

In Fig. 7, we present the energy loss function after the 
present BSE calculation and experiment in a wide energy 
range. Also for the onset region and the main three struc- 
tures below 20 eV the agreement is very close in positions 
and, a little less, in the oscillator strengths. For energies larger 
than 20 eV, a slight blue-shift of the peak positions is obtained 
with respect to the experimental ones. Furthermore, some dif- 
ferences with respect to the energy loss function —Im £ of 
other semiconductor/insulator systems could be considered. 
Indeed, these systems generally show a broad single structure 
around the plasmon frequency co p (e.g., as in the cases of cu- 
bic GaAs and Si),' 1 " unlike our theoretical curve (as shown 
in Fig. 7), while a structure of well-separated peaks appears 
in the spectrum of CdF^. A single broad structure would ap- 
pear around the energy of 34.6 eV, if all the valence electrons 
had contributed to a single plasmonlike excitation. Instead of 
this hypothetic structure, only two sharp peaks (located at 7.9 
and 13.4 eV, respectively) could be ascribed to excitonic ef- 
fects. Analyzing the positions of the peaks above 15 eV, it 
seems that they occur around the energies for which the real 
part of the dielectric function (see Fig. 3) shows its minima; 
i.e., around 19 eV, 25 eV, 30 eV, and 37 eV. In particular, one 



could suggest that due to the distinct separation of the valence 
states, in the case of CdFi the plasmon excitations take place 
but divided in well-defined energy windows. One could argue 
that the electrons remain at disposal for collective plasmonlike 
excitations as far as the corresponding energies are larger than 
the relative binding energies. Nevertheless, this hypothesis 
is correct for the F(2p), Cd(4d), and F(2s) valence electrons 
but it is no more true for the Cd(4p) and Cd(4s) electrons. In 
fact, considering the contribution of the F(2p) electrons one 
obtains a plasmon frequency of 20.53 eV. The contribution of 
F(2p) plus Cd(4d) electrons gives a plasmon energy of 27.81 
eV. Finally, a plasmon energy of 30.23 eV is obtained includ- 
ing the F(2s) electrons. These values are not so far away from 
the energies of the first three above-mentioned peaks of loss 
function. On the other hand, the Cd(4-p) and Cd(As) electrons 
lead to the following values for collective plasmonlike exci- 
tations: 32.53 eV [Cd(4p included] and 34.57 eV [Cd(4p) and 
Cd{As) included]. These energies are much smaller than the 
corresponding binding energies of about -62 eV for Cd(4p) 
and -102 eV for Cd{As) electrons. Therefore, a plasmonlike 
mechanism, involving Cd(4p) and Cd(4s) electrons, could not 
explain the structures in the high-energy spectrum of Fig. 7. 
It is interesting to note that the electron energy loss function 
of the rock salt rs-CdO, similarly to the present case, shows 
no evidence of a single pronounced plasma resonance coming 
from s, p or d electrons. 

In Fig. 8, the absorption coefficient a is shown in compar- 
ison with available experimental data. Considering the first 
part of the a function, the main four sharp peaks (at 7.8 eV, 
13.2 eV, 15.2 eV, and 18.0 eV) compare well with experimen- 
tal data. These facts can be ascribed to the four structures of 
the imaginary part of the dielectric function, taking in consid- 
eration that the absorption coefficient is directly proportional 
to the product of Im e(co) and the frequency, divided by the 
refractive index which shows a decreasing behavior in this en- 
ergy range. The three broader structures in a at higher ener- 
gies, between 20 eV and 40 eV, could be also posed in corre- 
spondence with three structures present in this energy range in 
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FIG. 8: Absorption coefficient after present calculation vs 
experiment (Ref. 3). 



the Im e(co) curve. For these structures a slight overestimate 
of around 0.5-1.0 eV occurs with respect to the experimental 
counterpart. 



and indirect gap of 1 1.8 eV). The Bethe-Salpeter scheme used 
for the optical absorption confirms the existence of an exciton 
located about 1.0 eV below the QP gap. The BSE spectrum 
for the imaginary part of the dielectric function shows an ex- 
cellent accordance with experiment. 

With respect to one-particle excitations we used the GW 
perturbative scheme for the electronic self-energy and calcu- 
lated the energy bands and the DOS. The electronic density of 
states at the gap region for CdFi and the energy-band structure 
have been compared with existing data in literature. The role 
of many-body effects turns out to be important in these one- 
particle properties by the large opening of the energy gaps. We 
show, moreover, that for the optical properties, many -body ef- 
fects, treated within the BSE scheme, are fundamental to ob- 
tain a reasonable comparison with existing experimental spec- 
tra. BSE imaginary part spectrum shows large differences in 
the shape which could not be recovered by any one-particle 
(GW/HSE) modified scheme at least in this material. 

Finally, the existence of an exciton located 1 eV below the 
quasiparticle gap for this compound has been discussed. 
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